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JH ' When extrasolar planets are observed to transit their parent stars, we are granted unprece- 

dented access to their physical properties. It is only for transiting planets that we are permitted 
direct estimates of the planetary masses and radii, which provide the fundamental constraints on 
models of their physical structure. In particular, precise determination of the radius may indicate 
the presence (or absence) of a core of solid material, which in turn would speak to the canonical 
formation model of gas accretion onto a core of ice and rock embedded in a protoplanetary disk. 
Furthermore, the radii of planets in close proximity to their stars are affected by tidal effects 
^ ■ and the intense stellar radiation. As a result, some of these "hot Jupiters" are significantly larger 

' than Jupiter in radius. Precision follow-up studies of such objects (notably with the space-based 

, platforms of the Hubble and Spitzer Space Telescopes) have enabled direct observation of their 

CO ' transmission spectra and emitted radiation. These data provide the first observational constraints 

CO ' on atmospheric models of these extrasolar gas giants, and permit a direct comparison with 

, the gas giants of the Solar system. Despite significant observational challenges, numerous 

transit surveys and quick-look radial velocity surveys are active, and promise to deliver an 
ever-increasing number of these precious objects. The detection of transits of short-period 
. Neptune-sized objects, whose existence was recently uncovered by the radial-velocity surveys, 

is eagerly anticipated. Ultra-precise photometry enabled by upcoming space missions offers 
^ , the prospect of the first detection of an extrasolar Earth-like planet in the habitable zone of its 

5_j ' parent star, just in time for Protostars and Planets VI. 

to 

^ 1. OVERVIEW structure and atmospheres may be compared in detail to the 

. planets of the Solar system, and indeed October 2005 was 

r/^ , The month of October 2005, in which the fifth Protostars . 1.1 r i. • • u - u i. r i. - . 

5-H notable tor being the month in which the number or objects 

' and Planets meeting was held, marked two important events • r . j 1 „ 

. . 1 o ' f jjjg former category surpassed the latter 

in the brief history of the observational study of planets or- „ . „ ^, . ... , . ^ u r , j 

■' J f Q^J• review or this rapidly-evolving field or study pro- 

biting nearby. Sun-like stars. First, it was the ten-year an- . ^ t c o j u • 1 

° -' ■' ceeds as follows. In Section 2, we consider the physical 

niversary of the discovery of 5 1 Pegb (Mayor fl«(i Qmc/oz, , , u- . u • ■ c^u 

■' ■' b \ J ^ ^ Structure of these objects, beginning with a summary of the 
1995), whose small orbital separation implied that similar 

hot Jupiters could be found in orbits nearly co-planar to our 

line of sight, resulting in mutual eclipses of the planet and 

star. Second, October 2005 heralded the discovery of the 



o 



observations (Section 2.1) before turning to their impact on 
our theoretical understanding (Section 2.2). In Section 3, 
we consider the atmospheres of these planets, by first sum- 
marizing the challenges to modeling such systems (Sec- 
ninth such transiting planet (BoMcny ef a/., 2005a). This se- T ,N J U .1 • • »U J . J 

^ ' tion 3.1), and subsequently reviewing the detections and 

lect group of extrasolar planets has enormous influence on i- j r wo t t\ 

^ ^ upper limits, and the inferences they permit (Section 3.2). 

our overall understanding of these obiects: The 9 transiting , , ... r . . /o a\ 

■> o We end by considering the future prospects (Section 4) 

planets are the only ones for which we have accurate esti- „, . ,^ , 1,1. j.uoi 
^ ■' for learning about rocky planets beyond the Solar system 

mates of key physical parameters such as mass, radius, and, ^, . j . j u . • c i, 1- ■ . 

■> ^ ■> ^ through the detection and characterization of such objects 

by inference, composition. Furthermore, precise monitor- • . -^^ c 

^ ^ in transiting configurations, 

ing of these systems during primary and secondary eclipse 

has permitted the direct study of their atmospheres. As a 
result, transiting planets are the only ones whose physical 
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2. PHYSICAL STRUCTURE 
2.1. Observations 

2.1.1. Introduction. When a planet transits, we can ac- 
curately measure the orbital inclination, i, allowing us to 
evaluate the planetary mass Mpi directly from the minimum 
mass value Mpi sin i determined from radial- velocity obser- 
vations and an estimate of the stellar mass, M^,. The plane- 
tary radius, Rpu can be obtained by measuring the fraction 
of the parent star's light that is occulted, provided a reason- 
able estimate of the stellar radius, i?*, is available. With 
the mass and radius in hand, we can estimate such criti- 
cally interesting quantities as the average density and sur- 
face gravity. Hence, the information gleaned from the tran- 
siting planets allows us to attempt to unravel the structure 
and composition of the larger class of extrasolar planets, 
to understand formation and evolution processes (including 
orbital evolution), and to elucidate physical processes that 
may be important in planetary systems generically. Fig. 1 
shows the mass-radius relation for the 9 known transiting 
planets, with Jupiter and Saturn added for comparison. It 
is fortunate that the present small sample of objects spans a 
moderate range in mass and radius, and appears to contain 
both a preponderance of planets whose structure is fairly 
well described by theory, as well as a few oddities that chal- 
lenge our present knowledge. 

We begin by describing how the objects shown in Fig. 1 
were identified and characterized, and, along the way, we 
illuminate the limitations that these methods imply for our 
efforts to understand extrasolar planets as a class. By def- 
inition, transiting planets have their orbits oriented so that 
the Earth lies nearly in their orbital plane. This is an uncom- 
mon occurrence; assuming random orientation of planetary 
orbits, the probabiUty that a planet with orbital eccentricity, 
e, and longitude of periastron, w, produces transits visible 
from the Earth is given by 
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which is inversely proportion to a, the orbital semi-major 
axis. All known transiting planets have orbital eccentricities 
consistent with zero, for which the last factor in the above 
equation reduces to unity. 

The radii of Jovian planets are typically only about 10% 
of the stellar radii. The transits known to date result in a 
0.3 — 3% diminution of the stellar flux reaching the Earth. 
These transits last for 1.5 — 3.5 hours, and accurate ground- 
based characterizations of these events are challenging. The 
paucity and subtlety of the transits make it necessary to use 
great care to reduce the random errors and systematic bi- 
ases that plague the estimation of the planets' fundamental 
properties (Section 2.1.4). 

2.7.2. Methods of Detection. The presently-known tran- 
siting planets have all been detected by one of the two fol- 
lowing means, both foreseen by Struve (1952): (1) Photo- 
metric detection of transit-like events, with subsequent con- 



firmation of planetary status via radial-velocity measure- 
ments, and (2) radial-velocity detection of a planet with 
subsequent measurement of photometric transits. Radial 
velocity detection has the advantage that the planetary na- 
ture of the target object is generally unambiguous. Its dis- 
advantage is that it requires substantial observing time on 
large telescopes to identify each planetary system, and only 
then can the relatively cheap process of searching for pho- 
tometric transits begin. Direct photometric transit searches 
simultaneously monitor large numbers of stars in a given 
field of view, but suffer from a very high rate of astrophysi- 
cal false positives (Section 2.1.3). 

Successful photometric transit searches have so far 
adopted one of two basic strategies, using either moderate- 
sized or very small telescopes to search either fainter or 
brighter stars. Five transiting planets (OGLE-TR-lOb, 56b, 
111b, 113b, and 132b) have been detected by the Optical 
Gravitational Lensing Experiment (OGLE) survey (Udal- 
ski et al, 2002a, 2002b, 2002c, 2003, 2004), which uses a 
1.3 m telescope. The parent stars of these planets are faint 
(typically V = 16.5). The large-telescope follow-up obser- 
vations needed to verify their planetary status, to measure 
the stellar reflex velocities, and to estimate the planetary 
masses and radii have been conducted by several groups 
{Bouchy et al, 2004, 2005b; Dreizler et al, 2002; Konacki 
et al, 2003a, 2003b, 2004, 2005; Moutou et al, 2004; Pont 
et al, 2004; and Torres et al, 2004a, 2004b, 2005). 

The Trans- Atlantic Exoplanet Survey (TrES) employed 
a network of 3 automated small-aperture (10 cm), wide- 
field (6° X 6°) telescopes {Brown and Charbonneau, 2000; 
Dunham et al, 2004; Alonso, 2005) to detect the planet 
TrES-1 {Alonso et al, 2004; Sozzetti et al, 2004). Its parent 
star {V = 11.8) is significantly brighter that the OGLE sys- 
tems, but fainter than the transiting-planet systems detected 
by radial-velocity surveys (below). Because of this relative 
accessibihty, TrES-1 has also been the subject of intensive 
follow-up observations, as detailed later in this review. 

Numerous other photometric transit surveys are active 
at the current time. The BEST {Rauer et al, 2005), HAT 
{Bakos et al, 2004), KELT {Pepper et al, 2004), Super- 
WASP {Christian et al, 2005), Vulcan {Borucki et al, 
2001), and XO {McCullough et al, 2005) surveys, and 
the proposed PASS {Deeg et al, 2004) survey all adopt 
the small-aperture, wide-field approach, whereas the EX- 
PLORE {Mallen-Ornelas et al, 2003) project employs 
larger telescopes to examine fainter stars. The benefits of 
surveying stellar clusters {Janes, 1996; Pepper and Gaudi, 
2005) have motivated several surveys of such systems, in- 
cluding EXPLORE/OC {von Braun et al, 2005), PISCES 
{Mochejska et al, 2005, 2006), and STEPSS {Burke et al, 
2004; Marshall et al, 2005). An early, stunning null result 
was the HST survey of 34,000 stars in the globular cluster 
47 Tuc, which points to the interdependence of the forma- 
tion and migration of hot Jupiters on the local conditions, 
namely crowding, metallicity, and initial proximity to O 
and B stars {Gilliland et al, 2000). 

Finally, three transiting planets were first discovered by 
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Fig. 1. — Masses and radii for the 9 transiting planets, as well as Jupiter and Saturn. The data are tabulated in Table 1, and are gathered 
from Bakos et al, in preparation, Bouchy et al. (2004, 2005b), Brown et ah, in preparation, Charbonneau et al. (2006), Holman et al. 
(2005), Knutson et al. (2006), Laughlin et al. (2005a), Moutou et al. (2004), Pont et al. (2004), Sato et al. (2006), Sozzetti et al. (2004), 
Torres et al. (2004a), and Winn et al. (2005). 



radial-velocity surveys. These include HD 209458b, the 
first transiting planet discovered (Charbonneau et al, 2000; 
Henry et al, 2000; Mazeh et al, 2001), and the two most 
recently discovered transiting planets, HD 149026b {Sato 
et al, 2005) and HD 189733b (Bouchy et al, 2005a). The 
latter two objects were uncovered by quick-look radial- 
velocity surveys targeted at identifying short-period planets 
of metal-rich stars (respectively, the N2K Survey, Fischer 
et al, 2005; and the Elodie Metallicity-Biased Search, da 
Silva et al., 2006). Given the preference of radial- velocity 
surveys for bright stars, it is not surprising that all three 
systems are bright (7.6 < V < 8.2), making them natural 
targets for detailed follow-up observations. As we shall see 
below, HD 209458b has been extensively studied in this 
fashion. Similar attention has not yet been lavished on the 
other two, but only because of their very recent discovery. 

2.1.3. Biases and False Alarms. Photometric transit 
surveys increase their odds of success by simultaneously 
observing as many stars as possible. Hence, their target 
starfields are moderately to extremely crowded, and the sur- 
veys must therefore work near the boundary of technical 
feasibility. The constraints imposed by the search method 
influence which kinds of planets are detected. 

Photometric transit searches are strongly biased in fa- 
vor of planets in small orbits, since such objects have a 
greater probability of presenting an echpsing configuration 
(Section 2.1.1). Moreover, most transit searches require a 
minimum of 2 (and usually 3) distinct eclipses to be ob- 



served, both to confirm the reality of the signal, and to per- 
mit an evaluation of the orbital period. Since larger orbits 
imply longer orbital periods and fewer chances for tran- 
sits to occur, small orbits are preferred for transit surveys 
with only a limited baseline. This is frequently the regime 
in which single-site surveys operate. However, multi-site 
surveys that monitor a given field for several months (e.g. 
HAT, TrES) frequently achieve a visibility (the fraction of 
systems of a given period for which the desired number of 
eclipse events would be observed) nearing 100% for peri- 
ods up to 6 days. As a result, such surveys do not suffer 
this particular bias, although admittedly only over a fim- 
ited range of periods. Similarly, a stroboscopic effect can 
afflict single-site surveys, favoring orbital periods near in- 
teger numbers of days and may account for the tendency of 
the longer-period transiting planet periods to clump near 3 
and 3.5 days (Pont et al, 2004, Gaudi et al, 2005). This 
situation occurs if the campaign is significantly shorter in 
duration than that required to achieve complete visibility 
across the desired range or orbital periods. However, for 
observing campaigns for which more than adequate phase 
coverage has been obtained, the opposite is true, and pe- 
riods near integer and half-integer values are disfavored. 
The limiting example of this situation would be a single- site 
campaign consisting of thousands of hours of observations, 
which nonetheless would be insensitive to systems with in- 
teger periods, if their eclipses always occur when the field 
is below the horizon. 

Most field surveys operate in a regime limited by the 
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signal-to-noise of their time series (which are typically 
searched by an algorithm than looks for statistically- 
significant, transit-like events, e.g. Kovdcs et ai, 2002), 
and for which the number of stars increases with decreas- 
ing flux (a volume effect). An important detection bias for 
surveys operating under such conditions has been discussed 
by Pepper et al. (2003) and described in detail by Gaudi et 
al. (2005), Gaudi (2005), and Pont et al. (2005). These 
surveys can more readily detect planets with shorter peri- 
ods and larger radii orbiting fainter stars, and since such 
stars correspond to a large distance (hence volume) they 
are much more numerous. As a result, any such survey 
will reflect this bias, which cannot be corrected merely by 
improving the cadence, baseline, or precision of the time 
series (although improving the latter will reduce the thresh- 
old of the smallest planets that may be detected). 

Most ongoing transit surveys are plagued by a high rate 
of candidate systems displaying tight curves that precisely 
mimic the desired signal, yet are not due to planetary tran- 
sits. We can divide such false positives into three broad 
categories: Some are true statistical false positives, result- 
ing from selecting an overly-permissive detection threshold 
whereby the light-curve search algorithm flags events that 
result purely from photometric noise outliers {Jenkins et al. , 

2002) . The second source is instrumental, due to erroneous 
photometry, often resulting from leakage of signal between 
the photometric apertures of nearby stars in a crowded field. 
However, the dominant form, which we shall term astro- 
physical false positives, result from eclipses among mem- 
bers of double- or multiple-star systems. Grazing eclipses 
in binary systems can result in transit-Uke signals with 
depths and durations that resemble planetary ones (Brown, 

2003) , and this effect is especially pronounced for candidate 
transits having depths greater than 1%. (For equal-sized 
components, roughly 20% of eclipsing systems have eclipse 
depths that are less than 2% of the total tight.) In these cases 
the ecUpse shapes are dissimilar (grazing eclipses produce 
V-shapes, while planetary transits have flat bottoms), but in 
noisy data, this difference can be difficult to detect. A false 
alarm may also occur when a small star transits a large one 
(e.g., an M-dwarf eclipsing a main-sequence F star). Since 
the lowest-mass stars have Jupiter radti, it is not surprising 
that such systems nunnic the desired signal closely: They 
produce flat-bottomed transits with the correct depths and 
durations. Larger stars ectipsing even larger primaries can 
also nninuc the desired signal, but a careful analysis of the 
transit shape can often reveal the true nature of the system 
(Seager and Mallen-Ornelas, 2003). Other useful diagnos- 
tics emerge from careful analysis of the tight curve outside 
of eclipses. These can reveal weak secondary eclipses, pe- 
riodic variations due to tidal distortion or gravity darken- 
ing of the brighter component, or significant color effects. 
Any of these variations provides evidence that the eclips- 
ing object has a stellar mass as opposed to a planetary mass 
(Drake, 2003; Sirko and Paczynski, 2003; Tingley, 2004). In 
the absence of these diagnostics, the stellar nature of most 
companions is easily revealed by low-precision (1 kms~^) 



radial velocity measurements, since even the lowest-mass 
stellar companions cause reflex orbital motions of tens of 
kms~^ (for examples, see Latham, 2003; Charbonneau et 
al, 2004; Bouchy et al, 2005b; Pont et al, 2005). 

The most troublesome systems are hierarchical ttiple 
stars in which the brightest star produces the bulk of the 
system's tight, and the two fainter ones form an eclipsing 
binary. In such cases, the depths of the eclipses are di- 
luted by tight from the brightest member, and often radial 
velocity observations detect only the bright component as 
well. Given neither radial velocity nor photometric evi- 
dence for a binary star, such cases can easily be mistaken 
for transiting planets. Correct identification then hinges on 
more subtle characteristics of the spectrum or light curve, 
such as tine profile shapes that vary with the orbital period 
(Mandushev et al, 2005; Torres et al, 2004b, 2005), or 
color dependence of the eclipse depth (O' Donovan et al, 
2006). Because of the large preponderance of false alarms 
over true planets, it is only after all of the above tests have 
been passed that it makes sense to carry out the resource- 
intensive high-precision radial-velocity observations that 
establish beyond question that the transiting object has a 
planetary mass. 

2.1.4. Determining the Radii and Masses. After transit- 
ing planets are identified, an arsenal of observing tools is 
available (and necessary) for their characterization. An ac- 
curate estimate of Mpi requires precise radial- velocity mea- 
surements (from which the orbital elements P, e, and ru 
are also determined), as well as an estimate of M^,. The 
former are gathered with high-dispersion echelle spectro- 
graphs fed by large telescopes. For bright parent stars, pre- 
cision of a few ms~^ (compared to reflex orbital speeds of 
50 — 200 m s~^) can be obtained with convenient exposure 
times, so that uncertainties in the velocity measurements do 
not donunate the estimate of Mpi . In this regime, the great- 
est source of uncertainty is the value of M^, itself. Given 
the difficulty of estimating the ages of field stars, compar- 
ison with grids of stellar models (e.g. Girardi et al, 2002) 
suggests that mass estimates are likely to be in error by as 
much as 5 % . This uncertainty could be removed by measur- 
ing the orbital speed of the planet directly. Several efforts 
have sought to recover the reflected-light spectrum of the 
planet in a series of high-resolution stellar spectra spanning 
key phases of the orbital period, but have achieved only 
upper Umits (Charbonneau et al, 1999; Collier Cameron 
et al, 2002; Leigh et al, 2003a, 2003b). (These results 
also serve to constrain the wavelength-dependent planetary 
albedo, a topic to which we shall return in Section 3.2.2.) 
For faint parent stars, the radial-velocity estimates become 
more expensive and problematic, and contribute signifi- 
cantly to the final error budget for Mpi. Interestingly, 
the most intractable uncertainty concerning masses of non- 
transiting planets, namely the value of sin i, is exquisitely 
well-determined by fits to the transit light curve. 

Analysis of moderate-precision light curves (obtained 
with ground-based telescopes) nonetheless yield a tight 
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Fig. 2. — //5r photometric light curves of transits of TrES-1 (top; Brown et al., in preparation) and HD 209458 (bottom; Brown et al. 
2001, offset by —0.007 for clarity). The shorter orbital period and the smaller size of the TrES-1 star result in a transit that is shorter in 
duration than that of HD 209458. Similarly, the smaller star creates a deeper transit for TrES-1, despite the fact that HD 209458b is the 
larger planet; the planetary sizes also affect the duration of ingress and egress. The TrES-1 data reveal a "hump" centered at a time of 
—0.01 d. This is likely the result of the planet occulting a starspot (or complex of starspots) on the stellar surface. 



constraint on the ratio Rpi/Ri,. However, fits to sucli data 
exliibit a fundamental degeneracy amongst tlie parameters 
Rpi, Ri,, and i, whereby the planet and stellar radii may be 
reduced in proportion so as to preserve the transit depth, and 
the orbital inclination may be correspondingly increased so 
as to conserve the chord length across the star. The uncer- 
tainty in Rpi is typically dominated by such degeneracies. 
Determining the value of Rpi requires fitting eclipse curves 
(facilitated by the analytic formulae of Mandel and Agol, 
2002) subject to independent estimates of M*, Ri,, and the 
stellar limb-darkening coefficients. If sufficient photomet- 
ric precision can be achieved, the value of iJ* may be de- 
rived from the light curve itself. This results in a reduced 
uncertainty on the value of Rpi, due to its weaker depen- 
dence on M^, (ARpi/Rpi) ~ 0.3(AM*/M*); see Char- 
bonneau (2003). For illustrative examples of the degenera- 
cies that result from such fits, see Winn et al. (2005), Hol- 
man et al. (2005), and Charbonneau et al. (2006). 

HST has yielded spectacular transit light curves for two 
bright systems, HD 209458 (Brown et al, 2001) and TrES-1 
{Brown et al, in preparation), which are shown in Fig. 2. 
The typical precision of these lightcurves is 10~* per one- 
minute integration, sufficient to extract new information 
from relatively subtle properties of the light curve, such 
as the duration of the ingress and egress phases, and the 
curvature of the light curve near the transit center. In prac- 
tice, such data have permitted a simultaneous fit that yields 
estimates of Rpi, R^,, i, and the stellar limb-darkening co- 
efficients, thus reducing the number of assumed parameters 



to one: M^. Cody and Sasselov (2002) point out that the 
combined constraint on {Mi^R^,) is nearly orthogonal to 
that resulting from light-curve fitting, serving to reduce the 
uncertainty in Rpi. Further improvements can result from 
the simultaneously fitting of multi-color photometry under 
assumed values for the stellar-Umb darkening, which serves 
to isolate the impact parameter (hence i) of the planet's path 
across the star and break the shared degeneracy amongst 
Rpi, R^, and i (Jha et al., 2000; Deeg et al., 2001). Re- 
cently, Knutson et al. (2006) have analyzed a spectrophoto- 
metric i/5rdataset spanning 290 — 1060 nm, and the com- 
bined effect of the constraints described above has been to 
permit the most precise determination of an exoplanet ra- 
dius to date (HD 209458b; Rpi = 1.320 ± 0.025 iJjup). 

2.1.5. Further Characterization Measurements. High- 
resolution stellar spectra obtained during transits can be 
used to determine the degree of aligiraient of the planet's 
orbital angular momentum vector with the stellar spin axis. 
As the planet passes in front of the star, it produces a char- 
acteristic time-dependent shift of the photospheric line pro- 
files that stems from occultation of part of the rotating stel- 
lar surface. This phenomenon is known as the Rossiter- 
McLaughlin effect (Rossiter, 1924; McLaughlin, 1924), 
and has long been observed in the spectra of eclipsing bi- 
nary stars. Queloz et al. (2000) and Bundy and Marcy 
(2000) detected this effect during transits of HD 209458. 
A fuU analytic treatment of the phenomenon in the context 
of transiting extrasolar planets has been given by Ohta et 
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al. (2005). Winnetal. (2005) analyzed the extensive radial- 
velocity dataset of HD 209458, including 19 measurements 
taken during transit. They found that the measurements 
of the radial velocity of HD 209458 during eclipse exhibit 
an effective half-amplitude of Av ~ 55 ms~^, indicat- 
ing a hne-of-sight rotation speed of the star of v sin = 
4.70 ± 0.16 kms~^. They also detected a small asymmetry 
in the Rossiter-McLaughlin anomaly, which they modeled 
as arising from an inclination, A, of the planetary orbit rel- 
ative to the apparent stellar equator of A = —4.4° ± 1.4°. 
Interestingly, this value is smaller than the A = 7° tilt of 
the solar rotation axis relative to the net angular momen- 
tum vector defined by the orbits of the solar system planets 
(see Beck and Giles, 2005). Wolfet al. (2006) carried out a 
similar analysis for HD 149026, and found A = 12° ± 14°. 
For these planets, the timescales for tidal coplanarization of 
the planetary orbits and stellar equators are expected to be 
of order 10^^ yr {Winn et al, 2005; Greenberg, 191 A; Hut, 
1980), indicating that the observed value of A likely reflects 
that at the end of the planet formation process. 

Pertubations in the timing of planetary transits may be 
used to infer the presence of satellites or additional plan- 
etary companions (Brown et al, 2001; Miralda-Escude, 
2002). Agol et al (2005) and Holman and Murray (2005) 
have shown how non-transiting terrestrial-mass planets 
could be detected through timing anomalies. Although HST 
observations have yielded the most precise timing measure- 
ments to date (with a typical precision of 10s; see tabulation 
for HD 209458 in Wittenmyer et al, 2005), the constraints 
from ground-based observations can nonetheless be used to 
place interesting limits on additional planets in the system, 
as was recently done for TrES-1 (Steffen and Agol, 2005). 

Precise photometry can also yield surprises, as in the 
"hump" seen in Fig. 2. This feature likely results from the 
planet crossing a large sunspot (or a complex of smaller 
ones), and thus is evidence for magnetic activity on the sur- 
face of the star. Such activity may prove to be an important 
noise source for timing measurements of the sort just de- 
scribed, but it is also an interesting object of study in its 
own right, allowing periodic monitoring of the steUar activ- 
ity along an isolated strip of steUar latitude (Silva, 2003). 

2.2. Theory and Interpretation 

2.2.7. Overview and Uncertainties. Transiting planets 
give us the opportunity to test our understanding of the 
physical structure of giant planets. In particular, structural 
models of the known transiting planets must be able to ac- 
count for the wide range of radiation fluxes to which these 
planets are subjected, and they must recover the observed 
range of radii. In general, as the planetary mass decreases, a 
given external energy input has an increasingly larger influ- 
ence on the size and interior structure of the planet. For hot 
Jupiters, the absorbed stellar flux creates a radiative zone in 
the subsurface regions that controls the planetary contrac- 
tion, and ultimately dictates the radius. Models of transiting 
giant planets straddle the physical characteristics of brown 



dwarfs and low-mass stars, as well as the solar system gi- 
ants (for an overall review, see Burrows et al, 2001). 

The construction of structural models for giant planets 
is difficult because a number of key physical inputs are 
poorly constrained. This situation holds equally for extra- 
solar planets and for the exquisitely observed outer planets 
of the Solar system. A benefit of robust determinations of 
the parameters for a growing range of planets is that un- 
certain aspects of the theory can become increasingly con- 
strained. Indeed, transit observations have the potential to 
clarify some of the core questions regarding giant planets. 

The dominant uncertainty regarding the overall structure 
of gas giants is in the equation of state (see the review of 
Guillot, 2005). The interiors of solar system and extrasolar 
giant planets consist of partially degenerate, partially ion- 
ized atomic-molecular fluids (Hubbard, 1968). The pres- 
sure, P, in the interiors of most giant planets exceeds 
10 Mbar, and central temperatures range from Tg ~ 10"* 
for Uranus and Neptune to Tc ~ 3 x 10^ for objects 
such as HD 209458b. This material regime lies beyond 
the point where hydrogen ionizes and becomes metalhc, 
although the details of the phase transition are still uncer- 
tain (Saumon et al, 2000; Saumon and Guillot, 2004). The 
equation of state of giant planet interiors is partially acces- 
sible to laboratory experiments, including gas-gun (Holmes 
et al, 1995), laser-induced shock compression (Collins et 
al, 1998), pulsed-power shock compression (Knudson et 
al, 2004), and convergent shock wave (Boriskov et al, 
2003) techniques. These experiments can achieve momen- 
tary pressures in excess of 1 Mbar, and they appear to be 
approaching the molecular to metallic hydrogen transition. 
Unfortunately, these experiments report diverging results. 
In particular, they yield a range of hydrogen compression 
factors relevant to planetary cores that differ by ~ 50%. 
Furthermore, the laboratory experiments are in only partial 
agreement with first-principles quantum mechanical cal- 
culations of the hydrogen equation of state (Militzer and 
Ceperley, 2001; Desjarlais, 2003; Bonev et al, 2004), and 
uncertainties associated with the equations of state of he- 
lium and heavier elements are even more severe (Guillot, 
2005). At present, therefore, structural models must adopt 
the pragmatic option of choosing a thermodynamically con- 
sistent equation of state that reproduces either the high- or 
low-compression results (Saumon and Guillot, 2004). 

Another uncertainty affecting the interior models is the 
existence and size of a radial region where hehum separates 
from hydrogen and forms downward-raining droplets. The 
possibility that giant planet interiors are helium- stratified 
has non-trivial consequences for their structures, and ulti- 
mately, their sizes. In the case of Saturn, the zone of helium 
rain-out may extend all the way to the center, possibly re- 
sulting in a distinct hehum shell lying on top of a heavier 
element core (Fortney and Hubbard, 2003). 

Uncertainties in the equation of state, the bulk compo- 
sition, and the degree of inhomogeneity allow for a de- 
pressingly wide range of models for the solar system giants 
that are consistent with the observed radii, surface tempera- 
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Table 1 

Properties of the Transiting Planets 



Planet 




P 






Rpi 


Mpi 


Teq pi 


Rpi 


Rpi 




Mq 


days 


K 






Mjup 


K 


20M© core 


no core 


OGLE-TR-56b 


1.04±.05 


1.21 


5970±150 


l.lOi.lO 


1.23±.16 


1.45±.23 


1800±130 


1.12±.02 


1.17±.02 


0GLE-TR-113b 


0.77±.06 


1.43 


4752±130 


0.76±.03 


1.081;°^ 


1.35±.22 


1186± 78 


1.07±.01 


1.12±.01 


UCjLIi- 1 K- 1 jzD 


1.35±.0o 


1.69 


6411±179 


1.43±.10 


1 1 o _i_ no 

1.13±.08 


1.19±.13 


1870±170 


1 1 O _l_ A'^ 

1.13±.02 


1 1 0_l_ A'^ 

1.18±.02 


HD 189733b 


0.82±.03 


2.22 


5050± 50 


0.76±.01 


1.15±.03 


1.15±.04 


1074± 58 


1.07±.01 


l.lli.Ol 


HD 149026b 


1.30±.10 


2.88 


6147± 50 


1.45±.10 


0.73±.05 


0.36±.03 


1533± 99 


0.98±.02 


1.15±.02 


TrES-1 


0.87±.03 


3.00 


5214± 23 


0.83±.03 


1.12±.04 


0.73±.04 


1038± 61 


1.02±.01 


l.lOi.OO 


OGLE-TR-lOb 


1.00±.05 


3.10 


6220±140 


1.18±.04 


1.16±.05 


0.54±.14 


1427± 88 


1.01±.02 


1.13±.01 


HD 209458b 


1.06±.13 


3.52 


6099± 23 


1.15±.05 


1.32±.03 


0.66±.06 


1314± 74 


1.02±.01 


1.12±.00 


OGLE-TR-lllb 




4.02 


5070±400 


0.85l:J° 


i.oot.JI 


0.53±.ll 


930±100 


0.97±.02 


1.09±.01 



tures, and gravitational moments. In particular (Saumon 
and Guillot, 2004), one can construct observationally- 
consistent models for Jupiter with core masses ranging from 
— 12 Mq, and an overall envelope heavy-element content 
ranging from 6 — 37 M^. This degeneracy must be broken 
in order to distinguish between the core accretion (Mizuno, 
1980; Pollack et ai, 1996; Hubickyj et aL, 2004) and grav- 
itational instability (Boss, 1997, 2000, 2004) hypotheses 
for planet formation. Fortunately, the growing dataset of 
observed masses and radii from the transiting extrasolar 
planets suggests a possible strategy for resolving the tangle 
of uncertainties. The extreme range of temperature condi- 
tions under which hot Jupiters exist, along with the variety 
of masses that are probed, can potentially provide definitive 
constraints on the interior structure of these objects. 

2.2.2. Comparison to Observations. Following the dis- 
covery of 51 Pegb (Mayor and Queloz, 1995), models of 
Jovian-mass planets subject to strong irradiation were com- 
puted (Lin, Bodenheimer and Richardson, 1996; Guillot et 
ai, 1996). These models predicted that short-period Jovian- 
mass planets with effective temperatures of roughly 1200 K 
would be significantly larger than Jupiter, and the discov- 
ery that HD 209458b has a large radius initially seemed to 
confirm these calculations. In general, Rpi is a weak func- 
tion of planet mass, reflecting the overall n = 1 polytropic 
character of giant planets (Burrows et ai, 1997, 2001). 

In order to evaluate the present situation, we have col- 
lected the relevant quantities for the 9 transiting planets in 
Table 1 . In particular, we list the most up-to-date estimates 
of P, Ri,, Mi,, Rpi, Mpi, as well as the stellar effective tem- 
perature, Tefi,-k- We also list the value of the planetary equi- 
hbrium temperature, Teq,pi, which is calculated by assum- 
ing the value for the Bond albedo. A, recently estimated for 
TrES-1 (A = 0.31 ± 0.14; Charbonneau et ah, 2005; Sec- 
tion 3.2.3). The precision of the estimates of the physical 
properties varies considerably from star to star. By drawing 
from the Gaussian distributions corresponding to the uncer- 
tainties in Table 1 and the quoted value for A, we can es- 
timate the uncertainty for Mpi and Teq,pi for each planet. 



Thereafter, for a particular choice of Mpi and Teq,p;, and 
fixing the planetary age at 4.5 Gyr, we can compute the- 
oretical radii. For this task, we use the results of Boden- 
heimer et al. (2003), who computed models for insolated 
planets ranging in mass from 0.11 — 3.0 Mjup. To evalu- 
ate the radii differences that arise from different heavy ele- 
ment fractions, separate sequences were computed for mod- 
els that contain and do not contain 20-M© solid cores, and 
both predictions are listed in Table 1 . The models have been 
calibrated so that, for the evolution of Jupiter up to the age 
of 4.5 Gyr, a model with a core gives the correct Jupiter 
radius to within 1%. Planetary age can also have a signif- 
icant effect on Rpi . For example, the evolutionary models 
of Burrows et al. (2004) for OGLE-TR-56b (1.45 Mjup) 
yield transit radii of Rpi ~ 1.5 i?jup at 100 Myr, and 
~ 1.25 i?jup after 2 Gyr. In general, however, Rpi evolves 
only modestly beyond the first 500 Myr, and hence the un- 
certainties in the ages of the parent stars (for which such 
young ages may generally be excluded) introduce errors of 
only a few percent into the values of Rpi . 

The models use a standard Rosseland mean photospheric 
boundary condition, and as such, are primarily intended for 
cross-comparison of radii. The obtained planetary radii are, 
however, in excellent agreement with baseline models ob- 
tained by groups employing detailed frequency-dependent 
atmospheres (e.g. Burrows et al, 2004; Chabrier et ah, 
2004; Fortney et al, 2005b). The models assume that 
the surface temperature is uniform all the way around the 
planet, even though the rotation of the planet is likely tidally 
locked. Hydrodynamic simulations of the atmosphere that 
aim, in part, to evaluate the efficiency with which the planet 
redistributes heat from the dayside to the nightside have 
been performed by Cho et al. (2003), Showman and Guil- 
lot (2002), Cooper and Showman (2005), and Burkert et al. 
(2005) under various simplifying assumptions. There is no 
agreement on what the temperature difference between the 
dayside and the nightside should be (Section 3.2.4), and it 
depends on the assumed opacity in the atmosphere. Burkert 
et al. (2005) suggest that with a reasonable opacity, the dif- 
ference could be 200 K, not enough to make an appreciable 
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difference in the radius. 

A number of interesting conclusions regarding the bulk 
structural properties of the transiting planets can be drawn 
from Table 1. First, the baseUne radius predictions display 
(la) agreement for seven of the nine known transiting plan- 
ets. Second, the planets whose radii are in good agreement 
with the models span the full range of masses and effective 
temperatures. The models do not appear to be systemati- 
cally wrong in some particular portion of parameter space. 
Although the reported accuracies of the basic physical pa- 
rameters are noticeably worse for the OGLE systems than 
for the brighter targets, the constraints are nonetheless use- 
ful to address models of their physical structure and, in par- 
ticular, the presence or absence of a solid core. Specifically, 
the baseline models in Table 1 indicate that the presence of a 
sohd core in a 0.5 Mjup planet with Tgs^pi = 1500 K leads 
to a radius reduction of roughly 0.1 i?jup. This difference 
generally exceeds the uncertainty in the estimate of Rpi . 

In the standard core-accretion paradigm for giant planet 
formation, as reviewed by Lissauer (1993), a Jovian planet 
arises from the coUisional agglomeration of a solid 10- M0 
core over a period of several million years, followed by a 
rapid accretion of hundreds of Earth masses of nebular gas, 
lasting roughly 10^ yr. The competing gravitational insta- 
bility hypothesis (e.g. Boss 1997, 2004) posits that gas-giant 
planets condense directly from spiral instabilities in proto- 
stellar disks on a dynamical timescale of less than 10^ yr. 
Boss (1998) points out that solid particles in the newly 
formed planet can precipitate to form a core during the ini- 
tial contraction phase. Only 1% of the matter in the planet is 
condensible, however, so a Jovian-mass planet that formed 
by this process will have a core that is much less massive 
than one that formed by the core-accretion scenario. 

Among the 7 planets that show agreement with the base- 
line models, it is presently difficult to discern the presence 
of a core. However, the "transit radius" effect {Burrows et 
al, 2003; Section 2.2.3) will tend to systematically increase 
the observed radii above the model radii listed in Table 1 
(which correspond to a 1-bar pressure level). Similarly, 
signal-to-noise-hmited field transit surveys bias the mean 
radius of planets so detected to a value larger than that of 
the intrinsic population (Gaudi, 2005). Taking both effects 
into account lends favor to the models with cores. Clearly 
more transiting planets and more precise determinations 
of their properties are necessary, as are more physically 
detailed models. We note that the identification of lower- 
mass transiting planets (for which the effect of a solid core 
is prominent) would be particularly helpful to progress in 
these questions. Several groups (e.g. Gould et al., 2003; 
Hartman et al., 2005; Pepper and Gaudi, 2006) have con- 
sidered the prospects for ground-based searches for planets 
with radii of that of Neptune, or less. 

2.2.3. The Transit Radius Effect. When a planet occults 
its parent star, the wavelength-dependent value of Rpi so 
inferred is not necessarily the canonical planetary radius at 
a pressure level of 1 bar (Lindal et al, 1981; Hubbard et 



al, 2001), which we have used for the baseUne predictions 
in Table 1 . As such, the measured radius is approximately 
the impact parameter of the transiting planet at which the 
optical depth to the stellar hght along a chord parallel to the 
star-planet line of centers is unity. This is not the optical 
depth in the radial direction, nor is it associated with the 
radius at the radiative-convective boundary. Hence, since 
the pressure level to which the transit beam is probing near 
the planet's terminator is close to one millibai (Fortney et 
al, 2003), there are typically 5 — 10 pressure scale heights 
difference between the measured value of Rpi and either 
the radiative-convective boundary (>1000 bars) and the 1- 
bar radius. (If, as discussed in Barman et al. (2002), the 
transit radius is at pressures well below the 1 millibar level, 
then the effect would be even larger.) Furthermore, exterior 
to the radiative-convective boundary, the entropy is an in- 
creasing function of radius. One consequence of this fact 
is significant radial inflation vis a vis a constant entropy 
atmosphere. Both of these effects result in an apparent in- 
crease of perhaps 0.1 i?,iup (^7%) in the theoretical radius 
for HD 209458b and 0.05 i?jup (~4%) for OGLE-TR-56b. 

2.2.4. Explaining the Oddballs. Two of the planets, 
HD 209458b and HD 149026b, have radii that do not agree 
at all with the predictions. HD 209458b is considerably 
larger than predicted, and HD 149026b is too small. These 
discrepancies indicate that the physical structures of the 
transiting planets can depend significantly on factors other 
than Mpi and Tca.pi . It would appear that hot Jupiters are 
imbued with individual personahties. 

While tiie radius of HD 209458b is certainly broadly 
consistent with a gas-giant composed primarily of hydro- 
gen, studies by Bodenheimer et al. (2001) and Guillot and 
Showman (2002) were the first to make it clear that a stan- 
dard model of a contracting, irradiated planet can recover 
Rpi ~ 1.35 iJjup for HD 209458b only if flie deep atmo- 
sphere is unrealistically hot. A number of resolutions to 
this conundrum have been suggested. Bodenheimer et al. 
(2001) argue that HD 209458b might be receiving interior 
tidal heating through ongoing orbital circularization. This 
hypothesis was refined by Bodenheimer et al. (2003), who 
computed grids of predicted planetary sizes under a vari- 
ety of conditions, and showed that the then-current radial 
velocity data set for HD 209458b was consistent with the 
presence of an undetected planet capable of providing the 
requisite eccentricity forcing. The tidal-heating hypothe- 
sis predicts that HD 209458b is caught up in an anomalous 
situation, and that the majority of hot Jupiter- type planets 
will have considerably smaller radii than that observed for 
HD 209458b. Recent analyses by Laughlin et al. (2005b) 
and Winn et al. (2005) indicate that the orbital eccentricity 
of HD 209458b is close to zero. This conclusion is further 
buttressed by the timing of the secondary eclipse by Deming 
et al. (2005a; discussed in greater detail in Section 3.2.3), 
which places stringent upper limits on the eccentricity, ex- 
cept in the unlikely event that the orbit is precisely aligned 
to our line of sight. Thus the eccentricity appears to be be- 
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Fig. 3. — Cut-away diagrams of Jupiter, Saturn, and the two oddball extrasolar planets, drawn to scale. The observed radius of 
HD 149026b implies a massive core of heavy elements that makes up perhaps 70% of the planetary mass. In contrast, the radius of 
HD 209458b intimates a coreless structural model, as well as an additional energy source to explain its large value. 



low the value required to generate sufficient tidal heating to 
explain the inflated radius. 

Guillot and Showman (2002) proposed an alternate hy- 
pothesis in which strong insolation-driven weather patterns 
on the planet drive the conversion of kinetic energy into 
thermal energy at pressures of tens of bars. They explored 
this idea by modifying their planet evolution code to include 
a radially adjustable internal energy source term. They 
found that if kinetic wind energy is being deposited at adi- 
abatic depths with an efficiency of 1%, then the large ob- 
served radius of the planet can be explained. Their hypothe- 
sis predicts that other transiting planets with similar masses 
and at similar irradiation levels should be similar in size to 
HD 209458b. The subsequent discovery that TrES-1 has 
a considerably smaller radius despite its similar tempera- 
ture, mass, and parent star metallicity is evidence against 
the kinetic heating hypothesis, since it is not clear why this 
mechanism should act upon only HD 209458b. 

Recently, an attractive mechanism for explaining the 
planet's large size has been advanced by Winn and Hol- 
man (2005) who suggest that the anomalous source of heat 
arises from obliquity tides that occur as a result of the planet 
being trapped in a Cassini state (e.g. Peale, 1969). In a 
Cassini state, a planet that is formed with a non-zero obliq- 
uity is driven during the course of spin synchronization to 
a final state in which spin precession resonates with or- 
bital precession. When caught in a Cassini state, the planet 
is forced to maintain a non-zero obliquity, and thus ex- 
periences continued tidal dissipation as a result of orbital 
libration. Order-of-magnitude estimates indicate that the 
amount of expected tidal dissipation could generate enough 
heat to inflate the planet to the observed size. 



HD 149026b presents a problem that is essentially the 
opposite to that of HD 209458b. Both the mass (0.36 Mjup) 
and the radius (0.73 i?jup) are considerably smaller than 
those of the other known transiting extrasolar planets. Cu- 
riously, HD 149026 is the only star of a transiting planet to 
have a metallicity that is significantly supersolar, [Fe/H] = 
0.36. The observed radius is 30% smaller than the value 
predicted by the baseline model with a core of 20 Mq. 
Clearly, a substantial enrichment in heavy elements above 
solar composition is required. The mean density of the 
planet, 1.17 gcni^"^, is 1.7x that of Saturn, which itself 
has roughly 25% heavy elements by mass. On the other 
hand, the planet is not composed entirely of water or sili- 
cates, or else the radius would be of order 0.4 or 0.28 i?jup, 
respectively (Guillot et ai, 1996, Guillot 2005). Models by 
Sato et al. (2005) and by Fortney et al. (2005b) agree that 
the observed radius can be recovered if the planet contains 
approximately 70 of heavy elements, either distributed 
throughout the interior or sequestered in a core. 

The presence of a major fraction of heavy elements in 
HD 149026b has a number of potentially interesting ramifi- 
cations for the theory of planet formation. Sato et al. (2005) 
argue that it would be difficult to form this giant planet by 
the gravitational instability mechanism (Boss, 2004). The 
large core also presents difficulties for conventional mod- 
els of core accretion. In the core-accretion theory, which 
was developed in the context of the minimum-mass solar 
nebula, it is difficult to prevent runaway gas accretion from 
occurring onto cores more massive than 30A/^, even if 
abundant in-falling planetesimals heat the envelope and de- 
lay the Kelvin-Helmholtz contraction that is required to let 
more gas into the planet's Hill sphere. The current structure 
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of HD 149026b suggests that it was formed in a gas-starved 
environment, yet presumably enough gas was present in the 
protoplanetary disk to drive migration from its probable for- 
mation region beyond 1 — 2 AU from the star inward to 
the current orbital separation of 0.043 AU. Alternately, a 
metal-rich disk would hkely be abundant in planetesimals, 
which may in turn have promoted the inward migration of 
the planet via planetesimal scattering {Murray et al., 1998). 

3. ATMOSPHERES 

By the standards of the Solar system, the atmospheres of 
the close-in planets listed in Table 1 are quite exotic. Lo- 
cated only 0.05 AU from their parent stars, these gas giants 
receive a stellar flux that is typically 10^ that which strikes 
Jupiter. As a result, a flurry of theoretical activity over the 
past decade has sought to predict (and, more recently, inter- 
pret) the emitted and reflected spectra of these objects (e.g. 
Seager and Sasselov, 1998; Seager et al., 2000; Barman et 
ah, 2001; Sudarsky et al, 2003; Allard et al., 2003; Bur- 
rows et al., 2004; Burrows, 2005). Observations promise 
to grant answers to central questions regarding the atmo- 
spheres of the planets, including the identity of their chem- 
ical constituents, the presence (or absence) of clouds, the 
fraction of incident radiation that is absorbed (and hence the 
energy budget of the atmosphere), and the ability of winds 
and weather patterns to redistribute heat from the dayside to 
the nightside. For a detailed review of the theory of extraso- 
lar planet atmospheres, see the chapter by Marley et al. We 
summarize the sahent issues below (Section 3.1), and then 
proceed to discuss the successful observational techniques, 
and resulting constraints to date (Section 3.2). 

3.1. Theory 

3.1.1. Overview. In order to model the atmospheres 
and spectra of extrasolar giant planets in general, and 
hot Jupiters in particular, one must assemble extensive 
databases of molecular and atomic opacities. The species 
of most relevance, and which provide diagnostic signa- 
tures, are H2O, CO, CH4, H2, Na, K, Fe, NH3, N2, and 
silicates. The chemical abundances of these and minority 
species are derived using thermochemical data and mini- 
mizing the global free energy. Non-equilibrium effects in 
the upper atmospheres require chemical networks and ki- 
netic coefficients. With the abundances and opacities, as 
well as models for the stellar spectrum, one can embark 
upon calculations of the atmospheric temperature, pressure, 
and composition profiles and of the emergent spectrum of 
an irradiated planet. With atmospheric temperatures in the 
1000 - 2000 K range, CO, not CH4, takes up much of 
the carbon in the low-pressure outer atmosphere, and N2, 
not NH;5, sequesters most of the nitrogen. However, H2O 
predominates in the atmospheres for both hot and cooler 
giants. Perhaps most striking in the spectrum of a close- 
in giant planet is the strong absorption due to the sodium 
and potassium resonance doublets. These lines are strongly 
pressure-broadened and likely dominate the visible spec- 



tral region. The major infrared spectral features are due to 
H2O, CO, CH4, and NH3. H2 collision-induced absorption 
contributes very broad features in the infrared. 

A self-consistent, physically reaUstic evolutionary calcu- 
lation of the radius, Tes^pi , and spectrum of a giant planet in 
isolation requires an outer boundary condition that connects 
radiative losses, gravity (g), and core entropy (S). When 
there is no irradiation, the effective temperature determines 
both the flux from the core and the entire object. A grid of 
TeS,pi > 9, and S, derived from detailed atmosphere calcu- 
lations, can then be used to evolve the planet (e.g. Burrows 
et at., 1997; Allard et al., 1997). However, when a giant 
planet is being irradiated by its star, this procedure must be 
modified to include the outer stellar fiux in the calculation 
that yields the corresponding S-Tes,pi -g relationship. This 
must be done for a given external stellar flux and spectrum, 
which in turn depends upon the stellar luminosity spectrum 
and the orbital distance of the giant planet. Therefore, one 
needs to calculate a new S-Tes.pi -g grid under the irra- 
diation regime of the hot Jupiter that is tailor-made for the 
luminosity and spectrum of its primary and orbital distance. 
With such a grid, the radius evolution of a hot Jupiter can 
be calculated, with its spectrum as a by-product. 

3.1.2. The Day-Night Effect and Weather A major issue 
is the day-night cooling difference. The gravity and interior 
entropy are the same for the day and the night sides. For a 
synchronously rotating hot Jupiter, the higher core entropies 
needed to explain a large measured radius imply higher in- 
ternal fluxes on a night side if the day and the night atmo- 
spheres are not coupled (e.g. Guillot and Showman, 2002). 
For strongly irradiated giant planets, there is a pronounced 
inflection and flattening in the temperature-pressure profile 
that is predominantly a result of the near balance at some 
depth between countervailing incident and internal fluxes. 
The day-side core flux is suppressed by this flattening of 
the temperature gradient and the thickening of the radia- 
tive zone due to irradiation. However, Showman and Guil- 
lot (2002), Menou et al. (2003), Cho et al. (2003), Burk- 
ert et al. (2005), and Cooper and Showman (2005) have 
recently demonstrated that strong atmospheric circulation 
currents that advect heat from the day to the night sides at a 
wide range of pressure levels are expected for close-in giant 
planets. Showman and Guillot (2002) estimate that below 
pressures of 1 bar the night-side cooling of the air can be 
quicker than the time it takes the winds to traverse the night 
side, but that at higher pressures the cooling timescale is far 
longer. Importantly, the radiative-convective boundary in 
a planet such as HD 209458b is very deep, at pressures of 
perhaps 1000 bar This may mean that due to the coupling 
of the day and the night sides via strong winds at depth, the 
temperature-pressure profiles at the convective boundary on 
both sides are similar This would imply that the core cool- 
ing rate is roughly the same in both hemispheres. Since the 
planet brightness inferred during secondary echpse (Sec- 
tion 3.2.3) depends upon the advection of stellar heat to the 
night side, such data can provide onstraints on the meteo- 
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rology and general circulation models. 

Almost complete redistribution of heat occurs in the case 
of Jupiter, where the interior flux is latitude- and longitude- 
independent. However, the similarity in Jupiter of the day 
and night temperature-pressure profiles and effective tem- 
peratures is a consequence not of the redistribution of heat 
by rotation or zonal winds, but of the penetration into the 
convective zone on the day side of the stellar irradiation 
(Hubbard, 1977). Core convection then redistributes the 
heat globally and accounts for the uniformity of the tem- 
perature over the entire surface. Therefore, whether direct 
heating of the convective zone by the stellar hght is respon- 
sible, as it is in our own Jovian planets, for the day-night 
smoothing can depend on the ability of the stellar insola- 
tion to penetrate below the radiative-convective boundary. 
This does not happen for a hot Jupiter. Clearly, a full three- 
dimensional study will be required to definitively resolve 
this thorny issue. 

Clouds high in the atmospheres of hot Jupiters with 
Tcs.pi > 1500 K would result in wavelength-dependent flux 
variations. Cloud opacity tends to block the flux windows 
between the molecular absorption features, thereby reduc- 
ing the flux peaks. Additionally, clouds reflect some of the 
stellar radiation, increasing the incident flux where the scat- 
tering opacity is high. This phenomenon tends to be more 
noticeable in the vicinity of the gaseous absorption troughs. 

3.2. Observations 

The direct study of extrasolar planets orbiting mature 
(Gyr) Sun-like stars may proceed without the need to im- 
age the planet (i.e. to spatially separate the light of the 
planet from that of the star). Indeed, this technical feat has 
not yet been accomplished. Rather, the ecUpsing geometry 
of transiting systems permits the spectrum of the planet and 
star to be disentangled through monitoring of the variation 
in the combined system light as a function of the known 
orbital phase. Detections and meaningful upper limits have 
been achieved using the following three techniques. 

3.2.1. Transmission Spectroscopy. The technique of 
transmission spectroscopy seeks to ratio stellar spectra 
gathered during transit with those taken just before or after 
this time, the latter providing a measurement of the spec- 
trum of the isolated star. Wavelength-dependent sources 
of opacity in the upper portions of the planetary atmo- 
sphere, or in its exosphere, will impose absorption features 
that could be revealed in this ratio. This technique can be 
viewed as probing the wavelength-dependent variations in 
the inferred value of Rpi. 

The first composition signature (Charbonneau et al, 
2002) was detected with the HST STIS spectrograph. The 
team measured an increase in the transit depth of (2.32 ± 
0.57) X 10"'* for HD 209458 in a narrow bandpass cen- 
tered on the sodium resonance lines near 589 nm. Ruling 
out alternate explanations of this diminution, they conclude 
that the effect results from absorption due to atomic sodium 



in the planetary atmosphere, which indeed had been unan- 
imously predicted to be a very prominent feature at visible 
wavelengths (Seager and Sasselov, 2000; Hubbard et al., 
2001; Brown, 2001). Interestingly, the detected amplitude 
was roughly 1/3 that predicted by baseline models that in- 
corporated a cloudless atmosphere and a solar abundance 
of sodium in atomic form. Deming et al. (2005b) follow the 
earlier work of Brown et al. (2002) to achieve strong up- 
per limits on the CO bandhead at 2.3 fim. Taken together, 
the reduced ampUtude of the sodium detection and the up- 
per limits on CO suggest the presence of clouds high in 
the planetary atmosphere (e.g. Fortney et al., 2003), which 
serve to truncate the effective size of the atmosphere viewed 
in transmission. Fortney (2005) considers the slant optical 
depth and shows that even a modest abundance of conden- 
sates or hazes can greatly reduce the size of absorption fea- 
tures measured by this technique. Alternately, non-LTE ef- 
fects may explain the weaker-than-expected sodium feature 
{Barman et al., 2005). 

Planetary exospheres are amenable to study by this 
method, as the increased cross-sectional area (compared 
to the atmospheres) implies a large potential signal. Vidal- 
Madjar et al. (2003) observed a 15 ± 4% transit depth of 
HD 209458 when measured at Lya. The imphed physical 
radius exceeds the Roche Umit, leading them to conclude 
that material is escaping the planet (Lecavelier des Etangs 
et al., 2004; Baraffe et al., 2004). However, the minimum 
escape rate required by the data is low enough to reduce the 
planetary mass by only 0.1% over the age of the system. 
More recently, Vidal-Madjar et al. (2004) have claimed 
detection of other elements, with a lower statistical signif- 
icance. Significant upper limits on various species at visi- 
ble wavelengths have been presented by Bundy and Marcy 
(2000), Moutou etal. (2001, 2003), Winn etal. (2004), and 
Narita etal. (2005). 

3.2.2 Reflected Light. Planets shine in reflected light 
with a visible-light flux fpi (relative to that of their stars, 

(?:),<«' ^(v)'«*^<«), 

where a is the orbital separation, p\ is the geometric albedo, 
and ^\{a) is the phase function, which describes the rel- 
ative flux at a phase angle a to that at opposition. Even 
assuming an optimistic values for px, hot Jupiters present 
a flux ratio of less than 10""' that of their stars. See Mar- 
ley et al. (1999), Seager et al. (2000) and Sudarsky et al. 
(2000) for theoretical predictions of the reflection spectra 
and phase functions of hot Jupiters. 

The first attempts to detect this modulation adopted a 
spectroscopic approach, whereby a series of spectra span- 
ning key portions of the orbital phase are searched for 
the presence of a copy of the stellar spectrum. For non- 
transiting systems, this method is complicated by the need 
to search over possible values of the unknown orbital inch- 
nation. The secondary spectrum should be very well sep- 
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Fig. 4. — The orbital phase-averaged planet-to-star flux-density ratio as a function of wavelength (A, in /nm) for the models of the 
four known transiting extrasolar planets for which such observations might be feasible (Burrows et al, 2005, and in preparation). 
The bandpass-integrated predicted values are shown as filled circles, with the bandwidths indicated by horizontal bars. The measured 
values for TrES-1 at 4.5 jim and 8.0 jj,m (Charbonneau et al, 2005) are shown as filled squares, and the observed value at 24 fim for 
HD 209458b {Denting et al., 2005a) is shown as a filled triangle. The extremely favorable contrast for HD 189733, and the extremely 
challenging contrast ratio for HD 149026, both result primarily from the respective planet-to-star surface area ratios. 



arated spectroscopically, as the orbital velocities for these 

hot Jupiters are typically 100 kms"'^, much greater than 
the typical stellar line widths of < 15 kms~^. Since the 
method requires multiple high signal-to-noise ratio, high- 
dispersion spectra, only the brightest systems have been ex- 
amined. A host of upper hmits have resulted for several 
systems (e.g. Charbonneau et al, 1999; Collier Cameron 
et al, 2002; Leigh et al, 2003a, 2003b), typically exclud- 
ing values of p\ > 0.25 averaged across visible wave- 
lengths. These upper limits assume a functional depen- 
dence for ^\{a) as well as a gray albedo, i.e. that the plan- 
etary spectrum is a reflected copy of the stellar spectrum. 

Space-based platforms afford the opportunity to study 
the albedo and phase function in a straightforward fashion 
by seeking the photometric modulation of the system light. 
The MOST satellite (Walker et al, 2003) should be able to 
detect the reflected light from several hot Jupiters (Green et 
al, 2003), or yield upper limits that will severely constrain 
the atmospheric models, and campaigns on several systems 
are completed or planned. The upcoming Kepler Mission 
(Borucki et al, 2003) will search for this effect, and should 
identify 100 — 760 non-transiting hot Jupiters with orbital 
periods of P < 7 d (Jenkins and Doyle, 2003). 

3.2.3. Infrared Emission. At infrared wavelengths, the 
secondary eclipse (i.e. the decrement in the system flux due 
to the passage of the planet behind the star) permits a de- 
termination of the planet- to- star brightness ratio. Since the 



underlying stellar spectrum may be reliably assumed from 
stellar models (e.g. Kurucz, 1992), such estimates afford the 
first direct constraints on the emitted spectra of planets or- 
biting other Sun-like stars. In the Rayleigh- Jeans Umit, the 
ratio of the planetary flux to that of the star is 

( IiL\ ~ -^eq,pi / Rpl \ ^ 

The last factor is simply the transit depth. From Table 1, we 
can see that the typical ratio of stellar to planetary temper- 
atures is 3.5 — 5.5, leading to predicted secondary eclipse 
amplitudes of several millimagnitudes. 

Charbonneau et al. (2005) and Deming et al. (2005a) 
have recently employed the remarkable sensitivity and sta- 
bility of the Spitzer Space Telescope to detect the ther- 
mal emission from TrES-1 (4.5 /xm and 8.0 pm) and 
HD 209458b (24 pm); Fig. 4. These measurements pro- 
vide estimates of the planetary brightness temperatures in 
these 3 bands, which in turn can be used to estimate (under 
several assumptions) the value of Tgq^pi and A of the plan- 
ets. Observations of these two objects in the other Spitzer 
bands shown in Fig. 4 (as well as the 16 pm photometric 
band of the IRS peak-up array) are feasible. Indeed, at 
the time of writing, partial datasets have been gathered for 
all four planets shown in Fig. 4. The results should per- 
mit a detailed search for the presence of spectroscopically- 
dominant molecules, notably, CII4, CO, and H2O. Using 
the related technique of occultation spectroscopy, Richard- 
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son et al. (2003a, 2003b) have analyzed a series of infrared 
spectra spanning a time before, during, and after secondary 
eclipse, and present useful upper limits on the presence of 
planetary features due to these molecules. 

Williams et al. (2006) have outlined a technique by 
which the spatial dependence of the planetary emission 
could be resolved in longitude through a careful monitor- 
ing of the secondary eclipse. Such observations, as well as 
attempts to measure the phase variation as the planet orbits 
the star (and hence presents a different face to the Earth) 
are eagerly anticipated to address numerous models of the 
dynamics and weather of these atmospheres (Section 3.1.2). 

The elapsed time between the primary and secondary 
eclipse affords a stringent upper limit on the quantity 
e cos w, and the relative durations of the two events con- 
strains esinro {Kallrath and Milone 1999; Charbonneau 
2003). The resulting limits on e are of great interest in 
gauging whether tidal circularization is a significant source 
of energy for the planet (Section 2.2.4). 

3.2.4. Inferences from the Infrared Detections. Varying 
planet mass, planet radius, and stellar mass within their er- 
ror bars alters the resulting predicted average planet-star 
flux ratios only slightly {Burrows et al., 2005). Similarly, 
and perhaps surprisingly, adding Fe and forsterite clouds 
does not shift the predictions in the Spitzer bands by an ap- 
preciable amount. Moreover, despite the more than a factor 
of two difference in the stellar flux at the planet, the predic- 
tions for the planet-star ratios for TrES-1 and HD209458b 
are not very different. Fortney et al. (2005b) explore the ef- 
fect of increasing the metallicity of the planets, and find a 
better agreement to the red 4.5/8.0 /imcolor of TrES-1 with 
a enrichment factor of 3 — 5. Seager et al. (2005) show 
that models with an increased carbon-to-oxygen abundance 
produce good fits to the HD 209458b data, but conclude that 
a wide range of models produce plausible fits. Barman et 
al. (2005) examine the effect of varying the efficiency for 
the redistribution of heat from the dayside to the nightside, 
and find evidence that models with significant redistribution 
(and hence more isotropic temperatures) are favored. 

In Fig. 4, there is a hint of the presence of H2O, since 
it is expected to suppress flux between 4—10 /um . This 
is shortward of the predicted 10 /um peak in planet- star flux 
ratio, which is due to water's relative abundance and the 
strength of its absorption bands in that wavelength range. 
Without H2O, the fluxes in the 3.6 — 8.0 /xm bands would 
be much greater. Hence, a comparison of the TrES-1 and 
HD 209458b data suggests, but does not prove, the pres- 
ence of water. Seeing (or excluding) the expected slope 
between the 5.8 /imand 8.0 /imbands and the rise from 
4.5 /umto 3.6 /um would be more revealing in this regard. 
Furthermore, the relative strength of the 24 /xm flux ratio in 
comparison with the 3.6 /im, 4.5 fim, and 5.8 /im channel 
ratios is another constraint on the models, as is the close- 
ness of the 8.0 (Umand 24 /xm ratios. If CH4 is present in 
abundance, then the 3.6 /imband will test this. However, 
the preliminary conclusion for these close-in Jupiters is that 



CH4 should not be in evidence. Models have difficulty fit- 
ting the precise depth of the 4.5 /im feature for TrES-1. It 
coincides with the strong CO absorption predicted to be a 
signature of hot Jupiter atmospheres. However, the depth of 
this feature is only a weak function of the CO abundance. 
A CO abundance 100 x larger than expected in chemical 
equilibrium lowers this flux ratio at 4.5 /im by only ~25%. 
Therefore, while the 4.5 /imdata point for TrES-1 implies 
that CO has been detected, the exact fit is problematic. 

In sum, Spitzer observations of the secondary eclipses of 
the close-in transiting giant planets will provide informa- 
tion on the presence of CO and H2O in their atmospheres, 
as well as on the role of clouds in modifying the planet- 
to-star flux ratios over the 3 — 25 fim spectral range. Fur- 
thermore, there is good reason to believe that the surface 
elemental abundances of extrasolar giant planets are not 
the same as the corresponding stellar elemental abundances, 
and Spitzer data across the available bandpasses will soon 
better constrain the atmospheric metallicities and C/O ratios 
of these planets. Moreover, and most importantly, the de- 
gree to which the heat deposited by the star on the day side 
is advected by winds and jet streams to the night side is un- 
known. If this transport is efficient, the day-side emissions 
probed during secondary eclipse will be lower than the case 
for inefficient transport. There is already indication in the 
data for HD 209458b and TrES-1 that such transport may 
be efficient (e.g. Barman et at., 2005), but much more data 
are needed to disentangle the effects of the day-night heat 
redistribution, metallicity, and clouds and to identify the 
diagnostic signatures of the chmate of these extrasolar gi- 
ant planets. The recently detected hot Jupiter, HD 189733b 
(Bouchy et al, 2005a) is a veritable goldmine for such ob- 
servations (Fig. 4), owing to the much greater planet- to- star 
contrast ratio. 

4. FUTURE PROSPECTS 

With the recent radial-velocity discoveries of planets 
with masses of 7 — 20 (e.g. Bonfils et al. 2005; Butler 
et al., 2004; McArthur et al, 2004; Rivera et al, 2005; San- 
tos et al, 2004), the identification of the first such object in 
a transiting configuration is eagerly awaited. The majority 
of these objects have been found in orbit around low-mass 
stars, likely reflecting the increased facility of their detec- 
tion for a fixed Doppler precision. Despite the smaller ex- 
pected planetary size, the technical challenge of measuring 
the transits will be alleviated by the smaller stellar radius, 
which will serve to make the transits deep (but less likely 
to occur). Due to the low planetary mass, the influence of 
a central core (Section 2.2.2) will be much more promi- 
nent. Furthermore, the reduced stellar size and brightness 
implies that atmospheric observations (Section 3.2) will be 
feasible. The radial- velocity surveys monitor few stars later 
than M4V, but transiting planets of even later spectral types 
could be identified by a dedicated photometric monitoring 
campaign of several thousand of the nearest targets. An 
Earth-sized planet orbiting a late M-dwarf with a week-long 
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period would lie within the habitable zone and, moreover, 
it would present the same infrared planet-to-star brightness 
ratio as that detected (Section 3.2.3). We note the urgency 
of locating such objects (should they exist), due to the lim- 
ited cryogenic hfetime of Spitzer. 

The excitement with which we anticipate the results 
from the Kepler (Borucki et al, 2003) and COROT (Baglin, 
2003) missions cannot be overstated. These projects aim to 
detect scores of rocky planets transiting Sun-like primaries, 
and the Kepler Mission in particular will be sensitive to 
year-long periods and hence true analogs of the Earth. Al- 
though direct follow-up of such systems (Section 3.2) with 
extant facilities appears precluded by signal-to-noise con- 
siderations, future facilities (notably the James Webb Space 
Telescope) may permit some initial successes. 

We conclude that the near-future prospects for studies of 
transiting planets are quite bright (although they may dim, 
periodically), and we anticipate that the current rapid pace 
of results will soon eclipse this review - just in time for 
Protostars and Planets VI. 
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